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Diffusion of a particle in the N- dimensional external poten- 
tial which is periodic in one dimension and unbounded in the 
other N — 1 dimensions is investigated. We find an analytical 
expression for the overdamped diffusion and study numeri- 
cally the cases of moderate and low damping. We show that 
in the underdamped limit, the multi-dimensional effects lead 
to reduction (comparing with the one-dimensional motion) 
of jump lengths between subsequent trapping of the atom in 
bottoms of the external periodic potential. As application 
we consider the diffusion of a dimer adsorbed on the crystal 
surface. 



I. INTRODUCTION 

A variety phenomena in physics and other fields can 
be modeled as Brownian motion in an external periodic 
potential One of particular example, the surface 

diffusion of atoms or small clusters, is of great fundamen- 
tal and technological interest During crystal growth 
the deposited atoms diffuse over the surface until they 
become incorporated in the lattice. On the semiconduc- 
tor Si(100) or Ge(100) surface, most of the deposited Si 
or Ge atoms combine to form dimers, and the diffusion of 
such dimers has recently been studied experimentally by 
a scanning tunneling microscope Moreover, atoms 
adsorbed on metal surfaces in some cases form closely 
packed islands which diffuse as a whole ||f7| . 

Theoretically the problem of diffusion can be described 
by a Langevin equation for the atom(s) or, equivalently, 
a Fokker-Planck-Kramers equation for the distribution 
function in the phase space In the trivial case of 

Brownian motion without the external potential, the dif- 
fusion coefficient is equal to D = Df = ksT '/m-q, where 
ks is the Boltzmann constant, T is the temperature, m is 
the particle mass, and r\ is the viscous damping coefficient 
which models the energy exchange with the substrate 
(thermostat). For a single atom in the one-dimensional 
(ID) sinusoidal external (substrate) potential the value 
of diffusion coefficient is well known as summarized in 
the Risken monograph Q. Exact results exist for the 
overdamped (Smoluchowski) case, r\ — > oo, when || 



D = D f l^ 2 (h) . 
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where h = e/2ksT, e is the (total) height of the substrate 
potential and lo(h) is the modified Bessel function, and 
in the underdamped limit, 77—5-0, when M 



where G(h) = (h/2ir) 1/2 e h J^ x (h) J(h), J{h) = 
f Q du u~ z l 2 e ~ 2h / u E _1 (u), and E(m) is the complete 
elliptic integral of second kind. At low temperatures, 
ksT <C e, both expressions (0) and (||) take the Ar- 
rhenius form, D = DA with A = cxp (— e/fc^T) and 
D u> 2 a 2 /2ttti in the high-friction case [here a is the pe- 
riod of the substrate potential and ujq — {2-rt / a)(e/2m) 1 / 2 
is the frequency of oscillation at its minimum], and 
D w irDf/2 in the low- friction limit. In a general case 
the diffusion coefficient can be found numerically with 
practically any desired accuracy by the matrix continued- 
fraction-expansion method Q . 

At low temperatures, ksT <C e, when the diffusion 
proceeds by uncorrelated thermally activated jumps over 
the barrier from one minimum of the external potential 
to another, the diffusion coefficient may be presented as 
D = RA(X 2 ), where RA is the rate of escape from a 
minimum of the external potential and (A 2 ) is the mean- 
square jump length. For a moderate or large damping, 
77 > wo, when only jumps for one period a of the exter- 
nal potential are possible, one should take A = a and 
R = RtstB(t]), where i?xST = ujo/2tt is the escape rate 
given by the transition state theory (TST) [|||]] , and the 
factor B(rj) = (z 2 + l) 1 ' 2 — z with z = rj/2uj s provides 
an interpolation between the TST and overdamped lim- 
its as was found by Kramers (T^| [here lo s is the "saddle" 
frequency at the saddle point x = x s , near which the ex- 
ternal potential has the form V(x) « e— |mUj(i — x s ) 2 ; 
for the sinusoidal potential uj s — luq] . The underdamped 
limit, 77 <^ ujq, is qualitatively different: in this case 
R w 2r\ej-nkBT oc 77 as was found firstly by Kramers JlO| , 
but the average jump length diverges as A oc rf~\ thus 
this again leads to the dependence i?A 2 oc 77" 1 similarly 
to the overdamped case. The occurrence of long jumps, 
A > a, has been observed in a number of experiments on 
surface diffusion jiJJT^]. The interval from low to moder- 
ate friction is covered by the Mel'nikov-Meshkov formula 
II 



Rmm ~ -Rtst exp I — 
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where A = Shri/too. Thus, the whole interval of fric- 
tions may be described by the interpolation formula 
R ~ -Rmm -8(77), which was checked numerically in [fl3|| . 
Combining this expression for R with the numerically 
calculated values of D, one can find the distribution of 
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jump lengths Q. Note that the widely used TST ex- 
pression D w i?TST^4a 2 , where the diffusion coefficient 
does not depend on the damping coefficient, operates in 
fact for a narrow interval of frictions close to the point 
rj ~ ujq only (which, fortunately, often corresponds to 
experimental situations). 

Although the described above results for one- 
dimensional diffusion are very important and often lead 
to reasonable estimations for experimentally measured 
diffusion coefficients, in real systems the motion always 
takes place in a N > 1 configurational space. Indeed, 
even for diffusion of a single atom adsorbed on a crystal 
surface N — 2 at least. Besides, the diffusing object may 
have internal degrees of freedom. Multi-dimensional ef- 
fects modify both the escape rate R and the jump length 
A. The escape rate can be presented as R — R\r>F, where 
the coefficient F is known as the "entropy factor" [fj"5| . 
The value of F can be found with the help of transi- 
tion state theory || which yields F w (ITwo,i) / 0^i^>s,i), 
where Wo,i are the frequencies at the minimum and u Si i 
are the "saddle" frequencies for all degrees of freedom i 
except the given one along the diffusion path. In this 
approach F can be interpreted as F = exp(AS l /fcs), 
where AS is the difference in entropy of the saddle and 
minimum-energy configurations. The entropy factor is 
often used to explain the "compensation effect" Q , when 
at experiment one observes that a decrease of the activa- 
tion energy (calculated as a slope of the Arrhenius plot 
of hiD versus T^ 1 ) is compensated by decreasing of the 
prefactor. As for the jump length, while for rj > u>q it still 
is given by A = a, in the underdamped limit it is mod- 
ified qualitatively comparing with the one-dimensional 
case. In the multi-dimensional space, the path connect- 
ing adjoining minima of the external potential may not 
coincide with the direction of easy crossing at the saddle 
point. Therefore, the probability of deactivation dur- 
ing long jumps is enhanced, leading to the reduction of 
jump length, A < Aid Jt6-18|. In particular, for the 2D- 
periodic substrate potential with the square symmetry it 
was found numerically [EsJ that D oc i]~ ' 5 which gives 
Acx?T - 75 . 

The multi-dimensional effects are also important in dif- 
fusion of molecules or small clusters: even for diffusion in 
the ID periodic potential (e.g., along "channels" on fur- 
rowed or stepped surfaces) one has for the dimer diffusion 
N = 2 at least. Diffusion of the dimer was studied numer- 



ically by Vollmer |19| with the help of matrix continued- 
fraction-expansion technique. The adiabatically slow mo- 
tion of a linear molecule in the ID sinusoidal potential 
was analyzed in ^o|, where the adiabatic trajectory was 
found for a general case. This allowed to find the activa- 
tion barriers and the minimum-energy and saddle-state 
frequencies and then to estimate the diffusion coefficient. 

The aim of the present paper is to study the multi- 
dimensional effects in diffusion processes. We consider 
two typical examples: motion of a single atom in a "chan- 
nel" which is periodic in one dimension and parabolic in 
others, and diffusion of a dimer (two-atomic molecule) 



in the ID sinusoidal potential. We find an analytical so- 
lution for the overdamped case and analyze numerically 
the dependence of diffusion coefficient on the damping 
constant rj. The numerical results were obtained with 
the Verlet algorithm by calculating the trajectory x(t) 
and then splitting it into N tT pieces, each of the time du- 
ration r. The diffusion coefficient was then calculated as 
D = (Ax 2 )/2t. 

The paper is organized as follows. In Sec. || we ob- 
tain the analytical expression for the diffusion coefficient 
in the overdamped limit. In Sec. Ill A we analyze the 
case of pure entropic barriers. In Sec. |[II B the activated 
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diffusion of a single atom is studied. In Sec._|II C the dif- 
fusion of a dimer is described. Finally, Sec. 
the paper. 



II. OVERDAMPED LIMIT 

Consider a particle moving in the A-dimensional ex- 
ternal potential Vn(x, yi, . . . , yN—i) which is periodic in 
the x direction, 



V N (x + a, . . .) = Vn(x, . . .), 



(4) 



and grows unboundently in the other N — 1 dimensions, 
V N (x, yi,..., y N - i) -> ]^ujf{x)yf if y l ->±co, (5) 

where u>i{x) > for all x. 

With presence of a viscous friction, the particle motion 
should be diffusive at long-time scale. The diffusion coef- 
ficient D can be found with the Einstein relation D = T/i, 
where the mobility // describes the proportionality be- 
tween the linear current j and the infinitesimal external 
dc force / which causes this current, j = jif . Therefore, 
we have to consider the particle motion in the external 
potential 



Vf(x,y u . . .,2/at-i) =V N - fx, 



(6) 



and then take the limit / — > 0. 

In the overdamped case, when the friction coefficient 
•q is much larger than the characteristic system frequen- 
cies, the motion of the particle is described by the Smolu- 
chowski equation 

dW - - - , 

— +V-J = 0, J = -r/~ (W Wf + TVW), (7) 

where W(x,yi, . . . ,yjv-i;i) is the distribution function, 
J(x,yi, . . . ,yjsr-i;t) is the density of particle's current, 
and the particle mass and Boltzmann constant are put 
to unity, m = 1 and ks = 1. 

For a steady state, Eq. ((?]) takes the form 



ox ox 



(8) 
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for the x component, and a similar form for other degrees 
of freedom. The density J of the current should satisfy 
the equation 



dJ x 
dx 



N a T 



(9) 



To reduce notations, below we consider the case of N = 
2 only; generalization to the N > 2 case is trivial. Let us 
introduce the one-dimensional density and current as 



p(x) 



4-oc 



dy W{x,y), 



+ OC 



dy J x (x,y). 



(10) 



(11) 



Owing to the condition the current j(x) does not 
depend on x, 



dj{x) 
dx 



-Jy(x, +00) + J V (X, — OO) = 0, 



(12) 



where we have used Eq. (Q). Thus, integrating both parts 
of Eq. (0) over y, we obtain the one-dimensional equation 



rp Mx) . / \ dV F {x) 
T — 1" P( x ) — — — = -VJ, 



dx 



dx 



(13) 



where we introduced the potential Vf(x) defined by the 
equation 



dV F (x) 



+ ~ 'dyW(x, y) d -Ylhyl. (14) 



dx J -co dx 

Now, if Vf (x) may be presented in the form 

V F (x) = V N (x;f)- fx, (15) 



where Vn(x; f) is a periodic function on x, Eq. (|13|) takes 
the form studied in and the diffusion coefficient can 
be calculated as 



D = £>/(/+/_)- 



I ± (T) = (2n)- 1 f W dxe ±v «^' T , 
Jo 

(16) 



where Df — T/rj and V c g(x) = lim/_>.o Vn(x; /). Thus, 
the diffusion coefficient D is determined by the one- 
dimensional function Vn(x;0). In the limit / — ► 
we may substitute the equilibrium distribution function 
W = W cq oc exp(-V N /T) into Eq. @, thus obtaining 



dV cS (x) $:™dye- v ^y)l T dV N (x,y)/dx 



dx 



JZ™ dy e- v »(*<y)/ T 



(17) 



Emphasize that this is the key approximation which 
is rigorous in the overdamped limit only. For the 



underdamped case, r\ — » 0, a similar multiplica- 
tive separation in the Fokker-Planck-Kramers equation, 
W(x,y,v x ,v y ,f) oc W(x,v x ,f)W e q(y,v y ), does not 
work even in the / — * limit. 
Let Vn(x, y) has the form 

V N (x,y) = V(x) + U(y)+v(x,y), (18) 

where the function v(x,y) describes the coupling between 
the two degrees of freedom. Then the effective potential 
V c g(x) can be presented in the following form, 



V eB (x) = V(x)-TS(x,T), 



(19) 



where the "entropy potential" S(x, T) is defined by the 
expression 

/ + CO 
dy exp{-[U(y)+v(x,y)}/T}. (20) 
-oo 

Notice that S(x) does not depend on V(x). 

III. APPLICATIONS 
A. Pure entropy barriers 



Let V(x) = in Eq. (18) 



and 



u (y) = 2 TOtJ iy 2 ' 



i 



v i x ,v) = T TO ( W 2 - w i) (1 - cosx)y , 



(21) 



(22) 



so that the atomic motion is inactivated in the x direc- 
tion, but the frequency of transverse oscillation depends 
on x, to = uj\ at x = and u> = L02 at x — tt. Then the 
integral in Eq. (^0|) can be easily evaluated analytically, 
and the entropy potential is given by the expression 



S(x) = - 1 -ln\l + 1 - 



^] -1 

UJi 



(1-cosx) } . (23) 



Notice that the entropy potential (g3|) does not depend 
on temperature, because both potentials ( f2l| ) and ( [22] ) 
depend on y in the same way (oc y 2 ). The function Six) is 
shown in Fig. |l]. It is periodic with the period a = 2n and 
the height es = | ln(w2/wi)|. The diffusion coefficient is 
given by D = DfF, where the entropy factor F depends 
on the ratio of frequencies z = ix>ijLO\ only, 



F(z) = [I + (z)I-(z)Y 



(24) 



where 



dx 



1 + -{z 2 - 1)(1 - cosa;) 



±1/2 



(25) 
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FIG. 1. The entropy potential S(x) for pure entropic bar- 
riers with wi/o»2 = 0.01 (dotted curve), toxjuj^ = 0.1 (solid 
curve), and lji/ui = 0.7 (dashed curve), respectively. 



reduces comparing with the free-diffusion value by a fac- 
tor of F (0.1) « 0.66 only. Indeed, although the height 
es tends to infinity at z — > 0, the width of barriers be- 
comes very narrow and thus cannot strongly modify the 
diffusion coefficient. 
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FIG. 2. The entropy factor F{w\/uto) for pure entropy 
barriers in the overdamped limit. 

Equ ations (HH yield F(z) = {tt/2) 2 K- 1 (VT^1 2 ) 
E _1 (-\/l — z 2 ), where K is the complete elliptic integral 
of first kind. Near z « 1 the function F(z) has the ex- 
pansion F(z) m 1 — |(1 — z) 2 , while at z — > it behaves 
as F(z) « (7r/2) 2 ln _1 (4/z). The function F(z) is pre- 
sented in Fig. 0. One can see that in the overdamped 
limit, the effect of entropy barriers is not too strong, in 
particular, even for ui/lu2 = 0.1 the diffusion coefficient 



FIG. 3. The diffusion coefficient D/D f (where 

Df — ksT/mri) as function of the friction coeffitient r\ for 
pure entropic barriers with uii = 0.1 and o>2 = 1 at T = 1. 
Inset: dependence on temperature (T = 1/3, 1 and 3) for 
r\ = 0.05 (up triangles), 0.5 (diamonds) and 5 (down trian- 
gles). 

In the underdamped case, on the contrary, the role of 
entropy barriers is essential. The dependence of the diffu- 
sion coefficient on the damping constant r\ was obtained 
numerically and shown in Fig. |3|. One can see that the 
function D(rj) exhibits a typical behavior of activated 
diffusion (D oc rj^ 1 at small and large frictions with a 
crossover between the limits) as might be expected from 
the shape of the entropy potential S(x) of Fig. |[ In 
the overdamped limit the average jump length is equal 
to the period of the potential S(x), X ~ 2ir, while in the 
underdamped limit long jumps with X/2tt 3> 1 play the 
dominant role as shown in Fig. ^ (in these simulations 
we assumed that the atom is trapped in a given well if 
it has sojourned in this well for a time lapse not shorter 
than (277) — 1 p|,pl|). The effect of entropy barriers is even 
stronger than might be expected from the analogy with 
the energy barriers of the same height. For example, for 
the frequencies uj\/lo2 = 0.1 used in the simulation, the 
height of the barrier is £5 = S(n) w 2.3, that would give 
the ratio D(t] -> oo)/D(r) -> 0) « 2e s /k B T » 4.6 for 
the T = 1 case, while the simulation leads to the ratio 
D(ri — > 00) /D(rj — > 0) > 33. From Fig. 0b one can see 
that (A/2tt) w 10 2 for the case of r] = 10 , while for the 
one-dimensional diffusion it should be (\/2ir) ~ r/^ 1 — 
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10 3 . Thus, multi-dimensional effects result in a strong re- 
duction of jump's length in the underdamped limit which 
leads to a decrease of the diffusion coefficient comparing 
with the fD motion. Note also that the dependence on 
temperature (shown in inset of Fig. ||) is almost negligible 
as has to be expected for the entropy potential. 



-Dsmoiuchowski = Dfl (e/2T). Because the entropy po- 
tential S(x) does not depend on the function V{x), it is 
still given by Eq. (E23h, and the integral ( |i"rj| ) can be eas- 
ily evaluated. The results for the overdamped limit are 
shown in Fig. ^, which can be compared with the simu- 
lation results for different frictions presented in Figures 
| and |. 
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FIG. 5. Activated diffusion with the barrier e = 2: D 
[normalized on the Smoluchowski value (|l|)] versus the ratio of 
transverse frequencies z = 0J2 /u>i in the overdamped limit for 
the temperatures T = 3 (dot-dashed curve), T = 1 (dashed 
curve), T = 1/3 (solid curve), and T = 1/9 (dotted curve). 
The short-dashed line shows the TST approximation. 



FIG. 4. Distribution of jump's length for pure entropic 
barriers with U\ = 0.1 and IV2 = 1 at T = 1 for (a) over- 
damped case (r/ = 1) and (b) underdamped limit (77 = 10 -3 ). 



B. Atom in a corrugated channel 

Let now the dependence of the external pot ential 
Vn(x,tj) on y is still given by Eqs. (|2l| ) and (p2[) , but 
the motion in the x direction is activated, 



1 



V(x) = -e(l 



cos a; 



(26) 



where e is the height of the external potential. At 
the minima of the potential ( p6| ) the transverse vibra- 
tions are characterized by the frequency oj\, while at 
the saddle points, by the frequency uj2- In the one- 
dimensional well as for the 2D case with lo\ = u>2, 
in the overdamped limit we have, according to Eq. (Q), 
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FIG. 6. D versus T for the activated motion with the 
barrier e = 2 for three values of the transverse frequencies 
(open diamonds and solid curves for u\ = u>2 = 1.0, down 
triangles and dotted curves for "wide barriers" with u>i — 1.0 
and ll>2 = 0.1, and up triangles and dashed curves for "narrow 
barriers" with lu\ — 0.1 and 012 = 1.0) for three values of the 
external damping: (a) r\ = 0.05, (b) r] — 0.5, and (c) 77 = 5. 
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From the D(T) dependence of Figures || and [6] one can 
see that at high temperatures, when the motion is in- 
activated, the uj\ = ll>2 case leads to the maximum of 
the diffusion coefficient similarly to the case with pure 
entropic barriers. With temperature decreasing, the en- 
ergy barriers and the entropy barriers play "in phase" 
for the "narrow-barriers" case of u)\ < u>2, and "in an- 
tiphase" for the "wide-barriers" case of uii > L02 ■ At low 
temperatures D > Did for the case of uj\ > u>2 at high 
and moderate frictions in agreement with predictions of 
the TST approach. The effect, however, is smaller than 
the TST predicts: in simulation we found that the dif- 
fusion coefficient changes only in three times when the 
ratio of frequencies is equal to ten. At very low frictions 
(e.g., r\ < 10~ 2 in Fig. ^), the entropy barriers become 
more important than the energy barriers, and the diffu- 
sion coefficient again becomes smaller than the ID one 
for all cases of lo\ ^ L02 as it was for the case of pure 
entropic barriers. 
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FIG. 7. The diffusion coefficient D [normalized on the 
Smoluchowski value ([!])] as function of the friction coefficient rj 
for the activated motion with the barrier e = 2atT = l/3 for 
three values of the transverse frequencies: (a) uji = L02 — 1.0 
[open diamonds, the dotted curves show the ID approximate 
values D « (ujo/2ir)Aa 2 B(ri) and D « ivDfA/2 at high and 
low frictions respectively], (b) cui = 1.0 and L02 = 0.1 (down 
triangles, "wide barriers"), and (c) ui\ = 0.1 and U12 = 1-0 (up 
triangles, "narrow barriers"). The dashed curves show the fit 

d{t)) oc ?r 1/3 - 

For moderate and low frictions the simulation results 
of Fig. can be fitted by a dependence D(rf) oc r/^ 1 ^ 3 . 
Because the escape rate R is still proportional to 77 in 
the multi-dimensional case p2|] , we may conclude that in 
the present case, the average jump length scales as (A) oc 
?7~ 2 / 3 , which is in agreement with the results of pure 
entropic barriers presented in the previous subsection, 



and also may be compared with the ID law (A) oc 77^ 1 
and the 2D simulation result (l8) (A) oc ?7~ 3 ' 4 . Thus, 
in the underdamped limit multi-dimensional effects lead 
to decreasing of diffusivity (comparing with the ID case) 
due to reduction of jump length which scales as (A) oc 
77 -2 / 3 instead of the ID scaling law (A) oc ry _1 . 



C. Diffusion of the dimer 

Now we can study diffusion of a dimer in the ID sinu- 
soidal potential. Let x\ and X2 arc the coordinates of two 
atoms coupled by the elastic spring with the constant g, 
and ao is the equilibrium distance (0 < ao < tt). Then 
the Hamiltonian of the system takes the form 

H = \^m a x\ + \n a x\ + ^ £ s(l - cos27ra;i/a s ) + 

ie s (f - cos27rx 2 /a s ) + i g (x 2 - x\ - a ) 2 . (27) 

In what follows we put e s — 2, m a = 1, a s — 2tt, and 
in the present paper we consider the case of ao = only. 
Introducing the coordinates x = x\ +X2 and y = X2 ~ x\, 
the Hamiltonian (|27]) can be rewritten as 



H= m (x 2 +y 2 ) + V N (x,y), 



(28) 



1 / x y \ 1 

V N {x,y) = -e (f - cos - cos -J + -gy , (29) 

which describes the motion of one particle of mass m = 
m a /2 = 1/2 in the a;-periodic potential of height e = 
2e a = 4 and period a = 2a s = Air. 

The adiabatic trajectory for this system was studied 
in pc| . Its shape depends on a value of the elastic con- 
stant g. The points (x, y) = (Aim, 0), where n is an inte- 
ger, always correspond to the absolute minimum of the 
potential energy. Near the minimum, the potential en- 
ergy has the expansion VN(x,y) « \m (u>q x x 2 + u>Q y y 2 ) 
with loqx = 1 and cjQ y — (2g + l) 1 / 2 . For a strong 
spring, g > 1/2, there is only one saddle point at 
{ x sjVs) = (27r, 0) between two adjacent minima (0,0) 
and (47r, 0). Near the saddle, the potential energy has 
the expansion 



Vn(x,v) 



1 



(30) 



with uj sx — 1 and uj sy = (2g — l) 1 / 2 , so that the acti- 
vation energy for dimer motion is equal to e s = e = 4 
(see Fig. ||). Therefore, dimer diffusion can be approx- 
imately described as motion of one atom in the corru- 
gated periodic potential with the transverse frequencies 
^1.2 = (2s±l) 1 / 2 , i.e. it corresponds to the case of "wide" 
barriers studied in the previous subsection. Thus, al- 
though the shape of adiabatic trajectory does not depend 



G 



on the elastic constant for the case of strong coupling, the 
diffusion coefficient does depend on g, it increases when 
g — » 1/2 due to decreasing of the transverse curvature at 
the saddle point. The simulation results of Fig. || show 
that the harmonic approximation describes the D(g) de- 
pendence with a good accuracy. From Fig. [h], where the 
ratio D(g)/D(Q) is presented for different temperatures, 
one can see also that close to the critical point g = 1/2, 
when anharmonicity of transverse vibrations at the sad- 
dle point is large, the entropy factor strongly depends on 
T, especially at low temperatures. 



(a) 




FIG. 8. (a) The activation energy e s and (b) the ratio of 
frequencies at the saddle and minimum points as functions of 
the elastic constant g for dimer's diffusion. 

For intermediate values of the elastic constant, 1/tt < 
g < 1/2, the adiabatic trajectory still has only one saddle 
point (27r, y s ) between the adjacent minima, where y s is 
now a solution of the transcendental equation sin(y s /2) = 
gy s . Near the saddle, the potential energy has the ex- 
pansion (Q) with the frequencies lu sx = [l — (gy s ) 2 ] ^ 4 
and uj sy — (2g — uj 2 ^ 1 / 2 . The saddle is characterized 
by the energy e s (g) = \e [1 + cos(y s /2)] + \gy 2 s , so that 
2 + tt/2 < e s < 4. 




U i I i I i I i I i u 

0.0 0.2 0.4 0.6 0.8 1.0 
g 

FIG. 9. The dependence of the diffusion coefficient D 
(times g) on the elastic constant g at T = 1 for different values 
of the damping constant: g — 5 (dotted diamonds), g = 0.5 
(open diamonds), g = 0.05 (solid diamonds), and g = 0.005 
(crossed diamonds). The dotted curve and plussed diamonds 
show the simulation results for the "atom in channel" model 
with g — 5 and other parameters adjusted to the dimer case. 

n 1 1 • 1 1 1 1 1 1 n 




g 

FIG. 10. The ration D(g)/D(0) as function of the elastic 
constant g for the dimer diffusion at g = 0.05 and different 
temperatures T = 3, 1, 1/2, and 1/3. 

Finally, for a weak coupling between dimer's atoms, 
g < 1/tt, there are two saddle points between the adja- 
cent minima (0,0) and (47r,0), with a local minimum of 
the potential energy between these saddle points. The 
coordinates of the saddle points are (27r — x',tt) and 
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(2tt + x' : 7r), where x' = 2 cos _1 (g7r). These saddle points 
are characterized by the energy e s (g) — ^(e + <?tt 2 ), 
so that 2 < e s < 2 + 7r/2. Near the saddle, the po- 
tential energy has the expansion (|3(]) with coefficients 
lu sx = (g - G) 1 / 2 and uo sy = (g + G) 1 / 2 , where G = 

[l-(- 2 '-D.9 2 ] 1/2 - 
The whole dependence s s (g) is shown in Fig. |8ja. 

The activation energy monotonically increases from the 
single-atom value e s = 2 at g = to the rigid-dimer value 
e s = 4 at g = 1/2 and then remains constant. Thus, one 
could expect that the diffusion coefficient should mono- 
tonically decrease with g increasing. However, the simu- 
lation results of Fig. [n] show that often this is not true. 
The peculiarity in the transverse frequencies at the point 
g = 1/2, where the saddle transverse frequency reaches 
zero, leads to a maximum of the function D(g) close to 
this point, if the damping is small, r\ < 0.5, and the tem- 
perature is not too low, T > 1 (recall e = 4). Thus, 
multi-dimensional effects may strongly affect dimer's dif- 
fusivity. 

IV. CONCLUSION 

In the present paper we studied in details the diffu- 
sion of a particle in two-dimensional space which is peri- 
odic along x and unbounded in the transverse direction. 
We calculated the entropy factor which emerges due to 
transverse degree of freedom, both in the overdamped 
limit (analytically) and in the undcrdamped case (nu- 
merically), and compared it with the prediction of the 
transition-state theory. We showed that in the under- 
damped limit, the multi-dimensional effects lead to re- 
duction (comparing with the one-dimensional motion) of 
jump lengths between subsequent trapping of the atom in 
bottoms of the external periodic potential. The simula- 
tion predicts that jump lengths scale as (A) / (Aid) <x r\ x l z . 
This leads to a decrease of diffusivity which now scales 
as D oc ry -1 / 3 instead of the ID dependence Did oc rf -1 . 

In the overdamped limit, the entropy factor (and, 
therefore, the prefactor in the Arrhcnius formula for acti- 
vated diffusion) does not depend on temperature as long 
as the transverse motion near the adiabatic trajectory 
may be described by the harmonic approximation. Sim- 
ulation shows that this remains true, at least approxi- 
mately, for low damping as well. Thus, in most cases 
the experimentally observed dependence of the prefactor 
on temperature has to be attributed to collective effects 
due to interaction between diffusing particles or/and be- 
tween the atom and (deformable) substrate. However, 
in the case of dimer diffusion at some value of the in- 
teraction between the atoms, when the saddle transverse 
frequency is equal zero, the anharmonicity of the trans- 
verse potential begins to play the important role and the 
entropy factor strongly depends on T. 
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